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Abstract 



Radiance is widely regarded as the principal quantity in light transport 
theory. Yet, the concept of radiance in use today has remained mostly 
unchanged since Lambert's work in the 18 th century. His formulation 
of the measurement of light intensity is based on classical differentials, 
and is known to suffer from several theoretical and practical limitations. 
After tracing the historic development of radiance and its shortcomings, 
we provide a modern formulation of light intensity measurements that 
models radiance as a differential form. We demonstrate the utility of 
this use of exterior calculus for questions in light transport theory by 
rigorously deriving the cosine term and the area formulation, without the 
need for postulates or heuristic arguments. The formulation of radiance 
as a differential form introduced in this paper hence provides the first 
step towards a modern theory of light transport. 



1 A Short History of Radiometry 

In ancient times, light intensity was believed to emanate from the eye as fire, one 
of the four elements. Already in the Middle Ages, however, it was realized that 
the intensity is an objective quantity that travels along light rays towards the 
eye. The first notable contributions toward a quantitative understanding of light 
intensity and its measurement were made in the 17 th century by Kepler and 
Mersenne who established the inverse square law and the cosine law for incident 
radiation. In the 18 th century, the work of BougueiQ and Lambert^] founded 
radiometry as a research field in its own right. Bouguer's work was concerned 
for example with the decay of light intensity in the atmosphere, which led him 
to what is now known as Beer-Lambert law. He also studied the scattering of 
light on rough surfaces modelled by specular microfacets. Lambert's central 
contributions in his Photometrica were: 

• the systematization of radiometry based on postulates firmly rooted in 
experiments; 
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• and the introduction of infinitesimals to quantitatively describe the mea- 
surement of light intensity. 

To model the measurement of light flux at a surface point in a direction, Lambert 
also introduced splendor, corresponding to luminance — the photometric sibling 
of radiance — in today's nomenclatur^] The physics of light played no role for 
Lambert, and in environments without absorption he employed an "actio in 
distans" principle to related outgoing to incoming intensity]^] Beginning in the 
18 th century, the relevance of radiometry for applications began to increase, and 
for example Euler employed it to study astronomical questions^] Lagrange in 
work on optical systems^] and Fourier in his seminal treatment of heat transfer^] 
The use of radiometry to study questions in opticsj^ astronomy^] and heat 
transport p*| continued throughout the 19 th century. Usually, however, the work 
employed only "the well-known laws of photometry"^] and even Lambert's 
Photometrica was rarely mentioned. Consequently, many important results 
were rediscovered over time. For example, the quadratic dependence of classical 
radiance on the refractive index was obtained independently for heatp^] and 
optical radiatiorp^] Although many prominent scientists and distinguished 
mathematicians, such as von Helmholtf^j Kirchhoff"^) and Clausiuf^j employed 
radiometry, the physical and mathematical foundations of the theory remained 
entrenched in the 18 th century. A notable exception is the work by LommeJ^j 
He considered radiance as emanating from a volume, an idea that was forgotten 
shortly afterwards, and also extended the Beer-Lambert law to include out- 
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scattering. Similar results that accounted for scattering were obtained by 
Chwolsorj^l and somewhat later by Schuster^] while Schwarzschilc p^l s tudied 
radiation equilibrium, a question again considered before by Lommap] Based 
on Schuster's workp^] Jacksorp] and Kin developed the transport equations 
that are in use to this date. The classical theory of radiometry culminated in 
the work by Nicodemusp^] In contrast to Lambert's postulates, he employed 
heuristic geometric arguments to establish the fundamental laws of radiometry. 
Nonetheless, conceptually Nicodemus' work belongs to the tradition of the 18 th 
and 19 th century. 

In the first half of the 20 th century, the introduction of methods from vector 
calculus led to a considerable improvement in the mathematical formulation 
of radiometric theor^J^] Vector irradiance led to closed form solutions for 
the irradiance for many geometric configuration^] and Arvo also employed 
the concept for light transport simulation^] As already recognised in 1939, a 
formulation of radiometry based on vector calculus is however limited since: 
"[vjaluable as the methods [using vector calculus] are, however, they probably 
do not constitute the ultimate solution of the problem. [. . . ] the physically 
important quantity is actually the illumination, which is a function of five 
independent, not three." p^] In the second half of the 20 th century, the physical 
foundations of radiometry and light transport were finally reconsidered. Early 
worl^Q] was based on particle models. This approach suffered however from an 
interpretation of light particles as photons which is physically not meaningful, 
cf. also the work by Mishchenko p*"] Later approached] were physically more 
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earnest. However, these efforts were limited by misconceptions about the nature 
of radiance, for example that it can only be expected to be meaningful in 
macroscopic environmental or that it represents an infinitesimal surface fluxp^] 
The 20 th century saw also many more applications of radiometry, for example 
in medical imaging, remote sensing, atmospheric science, astrophysics, and, last 
but not least, in computer graphics. However, the theoretical developments of 
the 20 th century were only rarely adopted by practitioners and most often the 
classical theory from Lambert's time was and is employed. 



2 Radiometry Reconsidered 

As we have seen, in the past 250 years radiometry has found applications in 
a wide range of fields — while most of the time it was not considered a subject 
worthy to be studied in its own right. The radiometric theory employed in 
practice is hence still that developed by Lambert in the 1750s. Central to this 
classical theory is the radiance 

L(x,uj) cos duo dA (1) 

which represents the infinitesimal flux of light energy in direction uj E S 2 
through a surface with normal n, with 6 being the angle between uj and np^] 
The importance of radiance lies in the ability to derive all other radiometric 
quantities — and the observable finite energy flux — from it by integration. De- 
spite this crucial role, however, the current conceptualization of radiance suffers 
from two central limitations: 

1. classical infinitesimals are employed; 

2. radiance is employed to model transport, although its mathematical 
formulation represents surface flux. 

As long as 80 years ago, these shortcomings led to a characterization of radiance 
as "absurdly antiquated"^] and in fact it is the current formulation that makes 
it necessary to introduce the cosine law and the area formulation as postulates or 
based on qualitative geometrical arguments, such as those commonly employed 
in the current computer graphics literature p] 

The current formulation of radiance using classical infinitesimals and its mis- 
use to describe transport has led to ample theoretical and practical limitations. 
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For example, for volume transport, where an "actio in distans" principle is not 
applicable, the length along the ray must be introduced a posteriori to obtain a 
volumetric description that can account for absorption or volume scattering. 
Likewise, in media with varying refractive index, an approach based on area 
fluxes is artificial. Because of this, the quadratic dependence of radiance on 
the refractive index can only be obtained using ad-hoc derivations such as inp^] 
Moreover, to describe transport in such inhomogeneous media, concepts from 
geometrical optics and radiometry have to be blended, even when no connection 
between them has been rigorously established to this datej^] Another significant 
problem with the current understanding of radiance is that differentiation along 
a ray is not possible. This prevents the use of derivative values in ray tracing, 
which would provide an avenue to better amortize the prohibitive costs of 
tracing rays by obtaining two or more values of the light energy signal from 
each single ray. For techniques that interpolate light energy, and which are 
frequently employed in practical applications, the intrinsic character of the light 
energy has to be respected: when interpolation is to be performed in space, 
such as in photon mapping, a volumetric formulation of radiance is required, 
while on surfaces, interpolation must ensure that measurements are correctly 
reproduced. 

In the present paper, we introduce a modern mathematical conceptualization 
of radiance that overcomes the limitations associated with the measurement 
of light intensity in classical radiometry. Our formulation of radiance as a 
differential form enables us to employ the geometric language of exterior calculus 
to rigorously establish the well-known cosine term and the area formulation, 
without the need for postulates or heuristic arguments as currently employed in 
the literature. Additionally, much insight into the nature of radiance is obtained. 
To avoid most of the complexity of modern light transport theoryp^] we restrict 
the following exposition to radiance and questions related to measurements. 
A comprehensive discussion, including a rigorous physical justification for the 
formulation of radiance we provide here, will be presented in a forthcoming 
publication. Some intuition is already provided in the appendix, and the 
intuition behind differential forms is provided in 



3 A Modern Formulation of Radiance 

Radiance describes the infinitesimal light energy flowing through a surface 
S C M 3 from direction p G S 2 in unit time and unit frequency. A finite 
measurement of energy flux is thus obtained by integration, which, to be 
physically meaningful, has to be covariant, e.g., independent of the coordinate 
system. Such covariant "integrands" are naturally modelled by differential 
forms-=a contemporary mathematization of classical infinitesimals. A modern 
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formulation of radiance can thus be written as 

A = L(q,p,v,t) dA^dpdv dt (2) 

where v and t are frequency and time, respectively, p = (pi,p2,P?>) is a unit 
direction (co-)vector and dp the infinitesimal solid angle form in this direction. 
The differential 2-form dA± e ft 2 (R 3 ) is given by 

dA± (p) = pi dy A dz + p2 dz A dx + ^3 A dy. (3) 

In contrast to the classical scalar- valued radiance, our formulation A of the 
concept is vector- valued, as is more apparent when written as 

A=( L(q,p, v)p\ dy A dz 

+ L(q,p,v) p2 dz A dx (4) 
+ L(q, p, v) ps dx A dy) dp dv 

or, when one omits the basis functions and just writes the components, as 

A = (L(q,p,i/)p u L(q,p,u)p 2 ,L(q,p,u)ps). (5) 

As will be seen in the next section, A becomes scalar only when pulled back to 
the chart of a surface, and it is this pullback that naturally leads to the cosine 
term. 

Remark 0.1. The most significant difference between our conceptualization 
of radiance and the classical one is the use of the 2-form dA± in Eq. [3| A 
2-form represents flux through an arbitrary surface as the linear superposition 
of fluxes through elementary 2-forms dy A dz, dz A dx, and dx A dy aligned with 
the coordinate axes. These elementary 2-forms are the 2-form basis vectors 
which represent the flux along the e x , e y , and e z axis, respectively, see the figure 
below. 
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Radiance differs from a classical, scalar 2- form in that it also has an angular 
dependence on p, and these two dimensions agree in that p represents the direc- 
tion of maximal flux, or the surface orientation for which the flux represented 
by dA± is maximized. 

Although vector- valued, our formulation of radiance should not be confused 
with vector irradiance J which is related to radiance through the integral 

J= / L(q,p,v) dpdv dA± dv dt, (6) 

where H 2 is the positive hemisphere over q G M 3 as defined by p. Eq. [6] implies 
that vector irradiance is also a differential 2-form, that is 

J = J(q)dA ± dvdt, (7) 

a fact which is also apparent when one re-reads Gershun's papei[^]with the 
modern theory in mind. The difference between vector irradiance and radiance 
is hence that the former is a "pure" 2-form, that is J G f£ 2 (R 3 ), whereas for the 
latter we have AgO 2 (R 3 )0 ^ 2 (T*Q), omitting in both cases the frequency 
and time dimension for clarity. 

In the following, we derive the cosine term and the area formulation directly 
from our formulation of radiance using the modern language of exterior calculus 
of differential forms. 



The Cosine Term 

A differential 2-form is an object that is naturally integrated over a surface 
5 C M 3 . As discussed for example in the supplementary material, this requires 
to pull back the form to a chart of the surface. For the radiance form A we 
thus have to determine 

[ A= [ far 1 )* A (8a) 

' A (^~'-Pl dudv ( 8b ) 
V (S) \du dv) 

where the local basis vectors d/du and d/dv for the tangent space TS are given 
by 

d/du = U x e x + U 2 e y + U 3 e z (9a) 
d/dv = V x e x + V 2 e y + V 3 e z (9b) 

and (U, (p) is a chart for 5, and we assume without loss of generality that all 
of S is covered by (U,cp), see Fig. [l] The pullback of A is most conveniently 
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Figure 1: Geometry of the pullback of the radiance 2-form A. 



computed using Eq. [8b] which, by linearity, can be determined term by term. 
With Eq. |4j we have for example for the first term 



(L(q,p, v)pi dy Adz) 



( d_ d_\ 
\du dv ) 



(10) 



where we also omitted the exterior basis functions unrelated to the surface 
integral. Again using linearity, it is sufficient to consider the basis function for 
the pairing and treat L(q,p,v)pi as a scalar factor. Eq. 10 can be simplified 
using 



(a A /3)(vi, v 2 ) = a(vi) f3(v 2 ) - a(v 2 ) fi(vi) 



(11) 



where a, (3 are arbitrary 1- forms and Vi,v 2 are arbitrary vectors]^] With Eq. [IT] 
we obtain 



(dy A dz) 



d_ d_ 
du dv 



--dy (U x e x + U 2 e y + U 3 e 3 ) dz (V x e x + V 2 e 2 + V 3 e 3 ) 
-dy (V x e x + V 2 e y + V 3 e z ) dz (U x e x + U 2 e y + U 3 e z ) 



and using the biorthogonality of the vector and covector basis functions this 
simplifies to 



(dy A dz) 



d_ d_ 
du dv 



U 2 dy(e y ) V 3 dz(e z ) - V 2 dy(e y ) U 3 dz(e z ) (12a) 
U 2 V 3 - U 3 V 2 . (12b) 
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Similarly, for the other two 2-form basis vectors we obtain 

( dzAdx )(^t) = uSvl ~ ulvS (13a) 

{dxAdy) (-L^) =uly2 - u2vl - (i3b) 

But the expressions in Eq. [l2]-Eq. |13b| are easily identified as the component 
form of the cross product. Since d/du and d/dv are the local basis vectors for 
the tangent space T<S, the pairing with the 2-form basis functions thus yields 
the components of the local surface normal, that is 



(14) 



With Eq. [TOj the integrand in Eq. [8b] can now be written as 

((/? -1 )*A = L(q,p, v) (pini J rp2n 2 +^3^3) dudv. (15) 

By writing the normal as n = ||n||n, where n = (ni, n2, ^3) is a unit vector, and 
with p ■ n = p\fii + P2U2 + P3TI3 we obtain 

((^~ 1 )*A = L(q,p, v)(p ■ n)\\n\\ dudv. (16) 
But by definition dA = \\n\\dudv so that 

(< f - 1 yA = L(q,p,v)(p-n)dA. (17) 

The vectors p and n are of unit length so that p ■ n = cos where is the angle 
between the vectors. We thus have for the infinitesimal area flux through the 
surface S in direction p that 

f d d\ 

A 7^—, dudv = L(q.p. v) cos dA dp dv. (18) 
\ou ov 



The cosine term in Eq. 18 appeared rigorously from our formulation of radiance 
as a differential 2-form without the need for postulates or heuristic arguments. 
It should also be noted that the dot product from which the cosine term was 
obtained was not derived from the Riemannian structure on M 3 but solely 
resulted from the pullback. 

Returning briefly to Eq. [l5j by linearity the equation can be written as 

L(q,p, v) (p • n) dudv dp dv = n • (L(g,p, v)pdpdv) dudv. (19) 

With the definition of vector irradiance in Eq. [6| we therefore have that the 
effective surface flux, the irradiance I(q), at a point q G M 3 is 



I(q)dt = n-\ / L(q,p,v)dp ) dA = n ■ J(q)dA (20) 





Figure 2: Geometry of the area formulation. 

where ^ C S 2 is now a finite solid angle, for example the solid angle subtended 
by a light source. Note that the possibility of computing the infinitesimal 
surface flux using a dot product was the original motivation for introducing 
vector irradianceF^I 



The Area Formulation of Light Transport 

In computer graphics, usually the parametrization of radiance over solid angle 
is employed. However, as was emphasized for example by Veachp^the area for- 
mulation is equally important]^] Historically it was in fact this parametrization 
that was used most frequently since it enables to obtain closed form solutions 
for the incident light energy. 

Central to the area formulation is the change of variable^] 

dp=^-dA(q) (21) 

relating the solid angle dp to an infinitesimal area dA(q) along the flow emanating 
from dp, see Fig. [2] To derive Eq. [21] we will employ the 2-form 

x dy A dz + y dz A dx + z dx A dy 
[x z + y z + z z ) 1 

where p = (x, y, z) and one should think of j3 G ^ 2 (M 3 \ {0}) as being centered 
at q G M 3 and "living" in the cotangent space T*R 3 ^ R 3 P] 

Remark 0.2. The 2-form ft G ft 2 (R 3 \{0}) is formed by 

1 



dA(p) (23) 
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where p has the natural interpretation as the direction of propagation. Equating 
p = n = (711,712,713), the surface 2-form dA G ^ 3 (M 3 \{0}) is given b$^\ 

dA(n) = -r-znr rii dy A dz + dz A dx -\- dx A dy. (24) 

HI 



The factor ofl/\\n\\ 2 in Eq. 23 thus ensures that the integral of ft over S 2 (r), the 
sphere with radius r, is always Att. This can be interpreted as always considering 
an equal size surface element of S 2 (r). 

We will begin by showing that j3 is a closed differential form, that is d/3 = 0. 
By linearity, we can compute the exterior derivative term by term, that is 

d/3=d ^w dy Adz ) +d \^ dz Ad ^ +d i^ dx A d y) (25) 

where ||p]| 3 = (x 2 + y 2 + z 2 ) 3 / 2 . For the first term of the exterior derivative we 
obtain 

x \ d x 
1, , lo J dy A dz = odx A dy A dz (26) 

\\p\\Zj dx( x 2 + y 2 + z 2f ^ > 

where the remaining partial derivatives vanish by the anti-symmetry of differ- 
ential forms. Using the product rule and after collecting terms, we then have 



2ij -\- -\- z 2 

dy Adz= — —-dx A dy A dz. (27a) 



JIpIIV (x 2 + y 2 + z 2 ) 5/2 

Analogously, for the other two terms of d/3 we obtain 



y \ y ^y 2 ~~\~ Z 2 

dz A dx — ^j-dx Ady Adz (27b) 



\\p\\ S J " ' (x 2 +y 2 + ^ 2 ) 5/2 

z \ y -\- y 2 — Iz 2 

dx Ady = ^Jodx A dy A dz. (27c) 



.\\P\\ 3 J " (x 2 +2/ 2 + z 2 ) 5/2 

With the above, it is immediately apparent that we have 

df3 = 5/2^ Ady Adz = (28) 

(x 2 + y 2 + z 2 ) ' 

and hence (3 G ^ 2 (R 3 \{0}) is indeed closed. 

In the following, we will also need the pullback of j3 G ^ 2 (IR 3 \{0}) onto an 
arbitrary surface S G M 3 . Following an argument analogous to those in Sec. [3J 
we obtain for the pullback 

{^YP=^dudv (29) 



49 Spivak, \Calculus on Manifolds: A Modern Approach to Classical Theorems of Advanced\ 
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where n is again the local surface normal. Writing p = \\p\\p and n = \\n\\n, 
where p and n are unit vectors, we obtain 



( '"»'' =H ir" = ir" (30) 



which is the right hand side of Eq. 21 Note that when the surface S coincides 
with a subset of the sphere S 2 so that the cosine term is unity and \\p\\ = 1, 
then we have = dp, that is, the pullback of p yields the solid angle 

measure. 

With d(3 = and Eq. 30 , we can finally establish the change of variables in 
Eq. 21 For this, consider a volume N C M 3 as shown in Fig. [3J bounded by 
the surface of a manifold M. and a subset of the sphere. Since /3 G f£ 2 (lR 3 \{0}) 
is closed, we have by Stokes theorem for differential forms that 

0= f d/3= [ (31a) 

JN JdN 

where d N is the boundary of N. Since dN consists of three parts, we can write 
= / /?+ / p + [ p. (31b) 

JdN 1 JdN 2 JdN 3 

The integral over dN2 vanishes since the radial vector p and the normal n of 
dN2 are orthogonal. We hence have 

= / p+ [ p. (31c) 

JdN 1 JdN 3 

By our previous results, the first term equals the solid angle subtended by M. 
while the second term is the integral over Eq. 30 Since M. was arbitrary, the 
result must also hold infinitesimally, that is, 

cos# 7 . . . 

dp = W ( } 



4 Discussion 

Radiance is the fundamental concept in radiometry, all other radiometric 
quantities being derived from it by integration. Nonetheless, its mathematical 
formulation still employs the classical infinitesimals that were introduced by 
Lambert. Moreover, while intrinsically related to measurements, it is today also 
employed to model the transport of light intensity. These inadequacies lead 
to numerous theoretical and practical problems, hampering the development 
of more effective computational techniques. In this paper, we provide a first 
step toward a modern theory of light transport by addressing the question of 
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Figure 3: Geometry employed in the proof of the area formulation of light 
transport. 



measurements and introducing a modern formulation of radiance. In contrast 
to classical radiance, our formulation as a differential form is vector-valued and 
becomes scalar only when paired to a surface, showing that radiance is the area 
flux through a yet undetermined surface. 

We demonstrated the utility of our geometric formulation of radiance by 
rigorously deriving the cosine term and the area formulation, results that previ- 
ously were either postulated or justified using heuristic geometric arguments. 
Our derivation also demonstrates that radiance falls off quadratically for mea- 
surements on concentric spheres when the area is kept constant, and that the 
classical results that "radiance is constant along a ray" has to be reconsidered, 
since radiance is not intrinsically a transported quantity, cf. the appendix. 

Readers familiar with vector calculus will have noticed that the result in 
Sec. [3] could have been obtained using classical means starting from a vector 
formulation of radiance, cf. Eq. [5] However, to the authors' knowledge, this is 
no onger the case for the derivation of the area formulation which requires a 
formulation of radiance as a differential form and the power of exterior calculus — 
also providing much added insight. Additionally, it is also the radiance 2-form, 
and not a vector, that can be derived from Maxwell's equations by considering 
measurements within the modern, geometric theory of light transport, as was 
developed in the first authors Ph.D. disseratiorf^l and will be presented in a 
forthcoming publication. 



Lessig, "Modern Foundations of Light Transport Simulation" 
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1 The Transport of a Density 



If radiance is only meaningful for measurements, how does one then describe the 
transport of light energy? To understand this, we will consider the transport 
of mass in fluid dynamics. Infinitesimally, the mass is described by the mass 
density p = a(q) dq G Den(Q), an element in Den(Q) which, roughly speaking, 
is equivalent to the space of volume forms Q 3 (Q) on Q C M 3 . The transport of 
the density in a fluid with velocity vector field u G X(Q) is then given by 

p + £np = 0, 

where £^p is the Lie derivative, or, in classical notation, 

p + V • (au) = 0. 




Let us now determine the mass flow through an infinitesimal surface sub- 
merged in the fluid. Using the transport theorem, a central result in continuum 
mechanic^] we obtain that the flow through a finite surface S C M 3 is 



(au) • n dA 



where ft is the surface normal. Infinitesimally, we thus have 

au - n dA 

which up to a factor that depends on the local magnitude of the fluid velocity 
vector field is equivalent to 

au • ndA — a cos dA, 

where u is the normalized velocity vector and = Z(u,n). When expressed 
using exterior calculus, the integrand pu is a 2-form, and given by p(q) dA u 
where dA u = u\ dq 2 A dq 3 + u 2 dq 3 A dq 1 + u$ dq 1 A dq 2 . 



51 Marsden, Ratiu, and Abraham, \Manifolds, Tensor Analysis, and Applications] Theorem 
.1.12. 
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As will be shown in a forthcoming publication, the quantity that is trans- 
ported in light transport theory is the phase space light energy density, or, 
short, light energy density, 

£ = L(q,p) dq dp G Den(T*Q) (33) 

which is a density on the six-dimensional phase space r*Q = QxR 3 forQcM 3 . 
The light energy density i G Den(T*Q) is the light transport equivalent of the 
mass density p G Den(Q), and the expression a(q) cos dA is the fluid equivalent 
of classical radiance L(x,lj) cos 6 duo d A, describing measurements through an 
infinitesimal surface. 
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